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Abstract 

We study the effective action for strong-coupling lattice QCD with one-component staggered 
fermions in the case of nonzero chemical potential and zero temperature. The structure of this 
action suggests that at large chemical potentials its ground state is a crystalline 'chiral density 
wave' that spontaneously breaks chiral symmetry and translation invariance. In mean-field theory, 
on the other hand, we find that this state is unstable. We show that lattice artifacts are partly 
responsible for this, and suggest that if this phase exists in QCD, then finding it in Monte-Carlo 
simulations would require simulating on relatively fine lattices. In particular, the baryon mass in 
lattice units, m^, should be considerably smaller than its strong-coupling limit of ~ 3. 

PACS numbers: 11.15.Ha,11.15.Me,25.75.Nq 
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I. INTRODUCTION 



The idea that the ground state of QCD at zero temperature and nonzero chemical poten- 
tial can spontaneously break translation invariance was first suggested by Deryagin, Grig- 
oriev and Rubakov (DRG) [1]. They studied large- QCD at weak coupling, and showed 
that the four-quark interaction, that is generated by a single gluon exchange, causes an 
instability in the quark Fermi sea toward the creation of a condensate of the form 

{ip{x)ip{x)) cos (2Q ■ x). (1.1) 

Here the wave vector Q has an arbitrary direction but its modulus is given by the Fermi 
momentum of the quarks 

Q=PF. (1.2) 

This phenomenon is often referred to as chiral density waves, and has famous solid-state 
counterparts like the one-dimensional Peierls instability [2j , and the spin density wave of the 
3- 



Overhauser effect 

The large- resu 
and Son showed in 



t of [li was consequently found to be misleading 0, Q], and Shuster 
^1 that if is decreased from infinity to < 0(1000), then color- 
superconductivity becomes preferable over the DGR state. Nonetheless, the authors in [6|] 
studied what may happen at lower densities (for A^ = 3), where weak coupling treatments are 
less reliable The conclusion of [6] was that chiral density waves may still be competitive 
with superconductivity, especially when one uses an instanton vertex to couple the quarks. 

Recently, there has been related progress in the Gross-Neveu model. This 1 + 1 model can 
be studied in large- A^, and its phase structure in the temperature-density plane was known 
to have a line of first order transitions that separate a low density phase with spontaneously 
broken chiral symmetry, and a high density phase, where this symmetry is intact. This 
structure arises when one restricts to an x-independent ansatz for {tp{x)ilj{x)) and was 
recently revisited in [8i] (and later with a lattice regularisation in [9|]), where the condensate 
was allowed to depend on the spatial coordinate of x. The study in discovered a new 
phase at low temperatures and high densities, which is separated by a second order phase 
transition from the low density phase, and where the chiral condensate has a crystalline 
structure. This crystal structure was also seen in the 't Hooft model p^], and is reminiscent 



of results from older work in the Skyrme model 11 1 . 
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In this paper we look for a phase with chiral density waves using an effective action that is 
derived from the strong-couphng expansion of lattice QCD. Apart from a phenomenological 
interest in this model, and in the way it may expose the phase we are after, we are also 
motivated by the following, more practical, reason. Future Monte-Carlo simulations of 
QCD at low temperatures and large chemical potentials, that will manage to control the sign 
problem, will presumably start by simulating relatively coarse lattices with strong couplings. 
This makes the knowledge on the phase diagram of the strong-coupling limit important, and 
in our context it is desirable to know whether one should expect a crystal phase, and if so, 
at which values of the chemical potential fi. 

Indeed, numerous authors have analytically investigated the phase diagram of the strong- 
coupling limit (for reviews see Q, Q), but the possibility to have chiral density waves was 
never considered, and these works were always restricted to homogeneous field configura- 
tions. 

Here we relax this constraint for the first time. To do so, we choose to work with one- 
component staggered fermions. In four dimensions, these lattice fermions describe four 
degenerate continuum Dirac fermions, commonly referred to as tastes. The continuum 
theory has an SU{A) x SU{A) chiral symmetry, that explicitly breaks on the lattice to the 
following axial taste non-singlet U{1) 

iIj ^ exp{i9-f5®^5) ip, ip ^ tp exp {i9 ^ ^5) . (1.3) 

Here ^5 is a 4 x 4 traceless matrix that operates in taste space, and that can be chosen 
to have the same matrix elements of 75. The ground state that we study in this paper is 
characterised by the helical condensates 

{tp{x) tp{x)) ~ cos (2Q ■ x) , (1.4) 
(^(x) (^ 75 ® ^5) ^(x)) ~ sin (2Q ■ x) , (1.5) 

that spontaneously break the chiral symmetry of Eq. (11.31) as well as translation invari- 
ance. Note that Eqs. f ll.4l) -( jL5l) cannot be rotated by a global transformation to have 
{ip{x) (75 (S> ^5) "ipix)) = 0. This means that if we regard the continuum limit as a theory of 



four degenerate light quarks, then our chiral waves break S'f/(4) flavor. As a result, we are 
studying a ground state which 
state, which is a flavor singlet. 



studying a ground state which is closer to the one studied in ^Ji] than to the original DGR 



The reason we choose this particular type of lattice fermions is that, in the strong-coupling 
limit, they give rise to a very simple effective action in terms of hadrons. This action was 
derived some time ago by Hoek, Kawamoto, and Smit [l^, and in the next section we re- 
derive it for self-completeness. After discussing why it is natural to expect chiral density 
waves in the ground state of this action (at nonzero /i), we formulate in Section UTTl a mean- 
field theory that is general enough to allow for such chiral waves to emerge. The mean-field 
equations are solved in Section IIVI where we also investigate how the mean-field ground 
state evolves when we increase fi. We conclude in Section |V] with a few remarks on the 
implications of this work and on future prospects. Appendix |A] includes technical details 
related to the calculation of the mean-field equations and the mean-field free energy. 



II. THE EFFECTIVE ACTION 



In this section we follow the work of Hoek, Kawamoto, and Smit (HKS) [15|], and re- 
derive their effective action. Readers who are familiar with [l^ can proceed to Section IIIBl 
where we discuss the structure of the effective action, and why it is natural to expect a chiral 
density wave in its ground state at nonzero fi. 



A. Deriving the Hoek-Kawamoto-Smit action for hadrons 

The starting point is the strong-coupling limit of an SU{N) lattice gauge theory with 
one-component staggered fermions. The action in this case is 

Xn] + E [NJnMn + C„5„ + 5„C„] , (2.1) 



71,1/ 



where n = (no, ni, . . . , n^) is an Euclidean lattice index and u = 0,1,2, ... ,d is the direction 
index [d is the number of spatial dimensions). The fields x and x are independent Grassmann 
variables and Unu is an SU (N) matrix that represents the gauge fields. The factors rjnu are 
the Kogut-Susskind factors that give the fermions their Dirac structure 

, jf" u = 0, 

Vnu = { (2.2) 

(_l)no+ni+...+n,_i ue[l,d]. 
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and here /i > is the quark chemical potential in units of the lattice spacing The 
source terms J, c, and c couple to the following meson and baryon fields 



1 ^ 

-^n = "77 Xn,aXn,ay 
a=l 

Bji = Xn,lXn,2 ■ ■ ■ Xn,N, 
Bn = Xn,NXn,N-l - ■ -Xn,!-, 



(2.3) 



and as indicated here, the only internal index that Xn and Xn carry, is the color index 
a = 1, . . . ,N. This leads to significant simplifications compared to the flavoured case, 
because it means that (M„)'' = for A; > A^, and (-Bn)'' = (Bn)'' = for A; > 1 30. 

Since the plaquette term is absent in Eq. fl2.ll) . one can readily integrate over the SU{N) 
link matrices. This is performed link by link and the result is 



Z(J,c,c) = jDxDxDUexpS = J DxDx exp S,{M, B, B) 



2-N+i 



+0) 



Bn ijlnu) Bn-\-u Bn-\-u 



N 



+ J2 {NJnMn + CnBn + 5nC„} 



Here the function Fjv is known for several values of N 18|, and in particular, for SU{3), it 
is 

(2.5) 



, . X 3 Q 15 q 

FqI-u) = -U-\ U U . 

' 4 64 256 



The next step is to write Z as a path integral over color-singlet fields. This is accomplished 
by first writing 



Z = exp 



Si{^dj,-dc,ds) 



Zq, 



Zo = J DxDxexp [nJ^M^ + B 
Performing the integral over x and x results in (here we assume that N is odd) 



(2.6) 
(2.7) 



^ + c c 



(2.8) 



that can be written as 



Zo = J DbDbDm exp |^ \^-dNbnm^^bn + NJ^nin + Cnbn + &nCn |, (2.9) 
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with dN = N\/N^. Here 6„ and 6„ are Grassmann variables, and the field m„ = e*^" takes 
values on the unit circle and has the measure 

/rf„„.i^=r* (2.10) 

We are now ready to apply Eq. fl2.6l) and doing so one obtains (setting J = c = c = 0, 
and replacing b ^ —2^^^h) 

Z = J DbDbDmexp{SnKs), (2.11) 

•S'hks = ^Meson + •S'Baryon + Sj, (2-12) 

•S'Meson = N Fn {mnnin+p) , (2-13) 
nu 

'S'saryon ^ ^ Dnm bmy (2'14) 
nm 

Si = Y.bnbn (2^'M^m-^). (2.15) 

n 

Here D is the massless Dirac operator of a free one-component staggered fermion whose 
chemical potential is equal to Nfi, 

d ^ 

^m,n+b^^^ ~ ^m,n-b^ + X! Vnu [^m,n+u " ^m,n-u] \ ■ (2.16) 



-Dn,m 2 



u=l 



Since the composite fields M„, i?„, and Bn in Eq. (12.11) couple to the same currents as m„, 
bn, and 6„ of S'hks (see Eq. (12.91) ). then their correlation functions are the same. This means 
that the fields 6„ and m„ represent the baryons and mesons respectively. Their hopping on 
the lattice is described by S'saryon and S'Meson, while the latter also describes meson-meson 
interactions. Interactions between mesons and baryons are described by the Yukawa-like 
vertex Sj, which in terms of quarks, represents a rearrangement of N quark-antiquark pairs 
into a single baryon-antibaryon pair. 



B. Chiral density waves from S'hks ? 

The purpose of this section it is to show why it is natural to expect that chiral density 
waves emerge from S'hks (Eq. (I2.12p ). We begin with Section [11 B II by formulating the chiral 
wave ansatz of Eqs. (I1.4l) - fll.5p in terms of the fields that appear in S'hks; and proceed to 
Section Tub 21 where we explain the mechanism by which these waves can arise. 
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1. The chiral density wave ansatz with staggered fermions 



We begin with constructing the condensates in Eqs. fll.4l) - fll.5l) . To do so we define a 
new lattice with spacing a = 2 that has 2'^^^ sites in its unit cell (This is the lattice of 
hyper-cubes used to define the taste basis pLS]). The original lattice coordinate n is related 
to the new one X by 

n = 2X + p, (2.17) 

where the {d + l)-dimensional vector p has py = 0, 1, and denotes the internal sites in the 
unit cell. The mean-field ansatz that we study in this paper is 

( (^)' ) = ( (rnnY ) = V, e'-^Q"-, g = 1, 2, . . . , iV. (2.18) 
where the sign factor e„ is given by 

Y^d „ +1 even site 
e„ = (-l)S-'"^= . (2.19) 

I —1 odd site 

The appearance of this sign factor in the phases of the condensates realises the helical 
structure of Eqs. fll.4l) - fll.5l) in the staggered formalism. To see this, we transform to taste 



basis with the unitary transformation U 19|] 



2^/2 2-D/2 

a=l f=l ^'^ '■" 

2D/2 2-D/2 

XpiX) = E E ^../(^)f^(a,/),p, (2.21) 

a=l f=l 

UiajiP = ^oIuyA . (2.22) 

\u=l ) a,f 

Here D = d + 1 is the number of spacetime dimensions, which we restrict to D = 1 + 1 or 
D = 3 + 1. The indices a and / are identified with the Dirac and taste indices respectively, 
and {a, /) is a composite index that takes 2^ values. The normalisation A/q is a chosen to 
have WU = 1, and the matrices 7^ are the Euclidean Dirac matrices.^ Using Eqs. fl2.20p - 
f l2.22p and the fact that the sign factor e„ depends only on p, 

en = e2x+p = {-l)^^"' =ep, (2.23) 



"'^ In d = 1 we choose (70,71,75) to be the Pauh matrices {(Jz,cfy,ax), and in c? = 3 we use the convention 
of [l9|. 
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we can write 



Y: {Xp{X)XpiX)) = Y: (i^aAXMX)) ^ (^(X)^(X)), (2.24) 



= {^^J{X)^ {^,®^,)tlj{X)). (2.25) 

Here 75 and ^5 have the same matrix elements, and act in Dirac and taste space respectively. 
Substituting Eq. fl238D in the left hand side of Eqs. ^Mj-^M) and using Eq. fl2:23|) gives 



{ipiX) iIj{X)) = A cos (2g ■ X + 0j , (2.26) 
(V'(X) {t 75 ® es) ^{X)) = -A sin (2Q ■ X + 0) , (2.27) 

with 

1 

= T^E^?- (2-28) 

{COS ihQl) d = 1 

^ , (2.29) 

cos(f) + T,l=i cos (^iQu - \T,,,^,yQu) d = 3 

By performing the chiral rotation of Eq. (11.31) with 6 = 0/2 we get rid of the angle in 

the arguments of the sine and cosine functions in Eqs. fl2.26p - fl2.27p . and obtain the assured 

helical structure of Eqs. fll.4p - fll.5l) . 



2. The mechanism of the chiral density wave instability 

We now proceed to show that the Yukawa-like interaction, Si, can make the state charac- 
terised by Eq. (12.181) become the ground state of S^ks (Eq. fl2.12p ) at /i > 0. For simplicity, 
we consider the d = 1 case only. 

In the absence of 5*/, the baryons behave like free massless fermions, and it is straight- 
forward to show that the poles of the propagator D^^ of Eq. (12.141) are determined in 
momentum space by 

= [sin^ (po/2 - ifiN) + sin^ (pi/2)]^ ; -tt < po,i < +7r. (2.30) 

Eq. fl2.30p tells us that ^Baryon describes the four energy bands 

sinhEj'^'^'^(pi) = +|sin(pi/2)|, (2.31) 
sinhEj^^'^^^(pi) = -|sin(pi/2)|, (2.32) 



with pi G [— TT, +7r], and that the Fermi energy and Fermi momentum of the baryons are 

Ep = fxN, (2.33) 
smpF/2 = sinhA^/i. (2.34) 

To take Si{b, b, m) into account we use mean-field theory, and replace the meson fields by 
their condensates. In the next section we find that with the mean-field ansatz of Eq. fl2.18p 
we should replace Si with 

^mean-field ^ ^ ^j^^ ^-iNQm ^ ^ ^2.35) 
n 

where the amplitude E is a complicated function of the amplitudes Vi^2,...,Af from Eq. (12.181) . 
In terms of the lattice coordinates of the new lattice Eq. (I2.35P becomes 

^mean-field _ ^5^(X)6^(X) (^Sc"'^ ^ (2^^+^^^ ^'') , (2.36) 
Xp 

where we defined bn = b2x+p = bp{X) and similarly for b. 

Clearly, 5'™ean-fieid j^j^gg baryons whose momenta differ by the amount 6pi = 2NQ. The 
strongest effect will occur between baryons that in the absence of Sj are degenerate in energy, 
and in particular, this mixing can occur between the bands E^'^^ and E^^) of Eq. flOTD . This 
will lead to level repulsion, and to the opening of a gap in the spectrum. To see this explicitly 
we move to momentum space with^ 

bp{X) = ./Ae-^^-^«(2Xi+pO/2 (2.37) 

V ^^s p 

bp(X) = ./Ae-p^Q(2^i+Pi)/2 ^e-*^^6p(p). (2.38) 

V p 

Note that Eqs. (I2.37p - p.38p are not a usual Fourier transform, because they gives the field 
bp{p) a momentum that differs from the momentum of bp{p) by an amount 6pi = 2NQep. 
This makes S'^'^^^~^'^^'^ diagonal in p space, and we find that the contributions of the baryons 
to ^KKS are given by 

^Baryon + ^r""""''" = HUp) Kpp'{p) bp\p), (2.39) 

p 

pp' 

^ Although we insert identical phases of epNQ{2Xi + pi)/2 into Eq. (|2.37p and Eq. (|2.38p . the transform 
between X-space and p-space still has a unit determinant due to the sign factor ep. (On a lattice with an 
even number of sites in each direction). 
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where the Dirac operator K is 



Kpp,{p) = [Fo sin (po/2 - iiiN) + Ti sin (pi/2 + eQN/2) + S 1]^^, . (2.40) 

Here 

(^),p' = (2.41) 

(r.)^^, ^ f,, (V,p+. + V,P-) e^^(^-''')/^ (2.42) 

with fjQ = 1 and r/i = (—1)'^". Using the methods of Appendix \K[ it is easy to calculate the 

determinant of K{p), extract its roots, and find that the four energy levels in Eqs. fl2.3ip - 
fl2.32p are modified to two positive energies, 



sinhiJ-i- = + 



S2 + (sin(pi/2) cos(iVg/2))^ ± cos(pi/2) sin(ArQ/2) 



(2.43) 



and two negative energies E± = —E±. We plot E± in Fig. [H where we choose S = Q = 
(black solid lines), S = 0.1 and NQ = 7r/4 (dashed blue lines), and E = 0.1 and (5 = 
(dotted green lines). Note that the case with S = corresponds to the absence of Si and the 
band structure should not depend on Q. Indeed, from Eq. (12.431) we see that if S = the 
energies E± reduce to (up to a shift of ±NQ from the origin, that can be removed 

by redefining the Brillouin Zone). 

Next, consider the following choice of the wave- vector Q 

NQ = pp. (2.44) 

Together with Eq. f l2.34p and Eq. f l2.3ip . this identifies the dashed-dotted horizontal line of 
Fig. [Hwith the Fermi level Eq. (12.331) . in which case the effect of the level repulsion is to 
lower all the energies in the Fermi sea (see dashed blue line in Fig. [1]). This decreases the 
baryon contribution to the free energy of the system with respect to the S = case. In 
contrast, if we choose Q = then Sj is a simple mass term that makes 



sinhE± = + sin2(pi/2), (2.45) 

and increases the energy of the Fermi sea (see dotted green line in Fig. [T]). As a result, if 
we denote the contribution of the baryons to the free energy by £^b(S, Q) then 

£B(S,p^/iV) < ^b(0,0) < ^b(S,0), (2.46) 
10 



-1 



-0.5 0.5 

-NQ/ti p/Jt +NQ/7I 



1 



FIG. 1: Dispersion relations E± : = Q = (black solid lines), S = 0.1 and NQ = 7r/4 (dashed 
blue lines), and T, = 0.1, and Q = (dotted green lines). Note that the solid and dotted lines have 
been shifted ±QN from the origin, so that they have the same abscissa as the dashed lines. With 
these shifts the level repulsion occurs at pi = 0, where the S = energy bands are degenerate with 
E+ = E_= sm{QN/2). 

and so, for /i > 0, the Fermi sea energetically prefers that the condensates in Eqs. f l2.18p 
carry a nonzero momentum Q. 

In the discussion above we have ignored many details. Firstly, we addressed only d = 1, 
where both the Fermi sphere and the gap occur at single points that can always be made 
to coincide. For d > 1, the Fermi sphere is a curved surface, and it is not assured that 
one can find a wave Q that will change the energy bands such that all the energies in the 
Fermi sea are pushed down. Indeed, it is known in condensed matter physics that similar 
inhomogeneous instabilities take place for d > 1 only for sufficiently strong interactions 

m 

. Secondly, we discussed only the contribution of the Fermi sea to the free energy. The 
remaining contributions are the anti-baryons energy (coming in the form of the negative 
energy bands E± of the Dirac sea) and the meson self-energy coming from S'Meson- In fact, 
as we find in the next sections, both these contributions depend on Q, and can prefer Q = 
for all values of /i > 0. Thirdly, S may also depend on Q, while in the discussion above we 
have treated it as being fixed. 
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Finally, note that our discussion relies on the level repulsion presented in Fig. [T], where 
we choose E = 0.1. Since Sj also gives mass to the baryons (see Eq. 02.451) ). this means 
that we have, so far, had in mind small baryon masses. It is a priori unclear whether the 
chiral wave instability occurs for large values of S and very massive baryons, but what is 
clear is that if the energy bands of the baryons in the crystalline phase look qualitatively 
like what we present in Fig. [H then it is natural to expect this instability when NQ ~ Pi?. 
More precisely, by comparing the solid black lines and the dashed blue lines in Fig. [H we 
see that what allows the crystalline phase (that has nonzero S and Q) to compete with the 
massless phase (that has S = 0) is the fact that it has gapless excitations. 

To study the effects of the issues we mention above on the viability of the chiral density 
waves, we develop a mean-field analysis for inhomogeneous vacua in the next Section. 

III. MEAN-FIELD THEORY 



We start from Eq. (12.121) and introduce auxiliary fields for the expectation values of 
rrin, (m„)^, . . . , (m„)^ by writing 



exp 



exp 



N 



q=l nu 



q=l nu 



KVq,n - ml) 



(3.1) 



= j DV^exp 
= j DVDhexp 

Here / dV = Y[ (^^r^ / dVj''\ and V^j are the real and imaginary parts of Vq, 

J q J — oo J — oo ' 

dh = Y[ dh^^^ / dh^^\ Also, hqVq means h^^V^^ + h'^^^vi'^\ and similarly 

„ J —ioo J —ioo 



^ E E «gK,nK,n+^ - Nhq,^ (Vq,^ - m^) 
5=1 nu 



for hq{mY. Proceeding from Eq. (13.11) . the action for the scalar fields V and h becomes 



N 



Sesiy, h) = E '^9 E ^q,nVq,n+u - N hq,nVq,n + Sq, 
q=l nu 



(3.2) 



n,q 



exp(5o) = / DbDbDmexpl b- h^-^dNm-^l + d] ■ b + N J2 Kn^l} ■ (3.3) 



q,n 



The starting point of Mean field theory is to write the path integral for the action in 
Eq. as 

Z{X) = J DhDVe^^^^ (3.4) 
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and to approach the physical point, A = 1, by expanding around A = [20|. Calculating 
0(A) corrections is, however, outside the scope of the work we present here, and we study 
only the 0(A°) level. In this case, the tree-level mean- field equations are 



-— = 0, and ——-^ = 0, (3.5) 



which give 



±d 

= E Vg,n+i>, (3.6) 

u=±0 
±d 

hfl = ^a, Y: V<i,n+o, (3.7) 

v=±0 

Vg,n = «). (3.8) 

Here the average (, ) in Eq. (13. 8p is with respect to the path integral in Eq. (13.31) . 

Equations fl3.6p - fl3.7l) tell us to try an ansatz with h^'^'^ = ih^^^ = ih, in which case the hV 
and hm terms in Eq. (13.21) and Eq. (13. 3p have the meaning of a simple complex multiplication. 
Restricting to this type of ansatze we use 

dm 



w!" = 4,0, (3.9) 



to perform the J Dm integral in Eq. (13.31) . and we find (for = 3) 

/ dm r / , 2 7 •hX 8 n- 

i> exp 3 him + h2m + h^m H — m bb 

J 2'Kim L ^ ^9 



exp 



A[h\ + 2hih2 + ^/i3 ) bb 



(3.10) 



This turns Z{1) of Eq. (EJ) to 

Z = J DV DhDbDb expSes, (3.11) 

N 

SeS = NY^agY. Vg,nVg,n+u ' + E 1 + D]^^ b^, (3.12) 
g=l ni/ n,q n,m 

Sn = 4 (^/l? „ + 2/li,„/l2,„ + ^/l3,n) • (3.13) 

The next step is to assume an ansatz for the spatial behaviour of V^_„ and hq^n- Here we 
allow for the possibility of the broken translation invariance, and write 

Vq^n = Vq CXp {+iq tn Q uj , /ig,n = kg CXp (-ig tn Q uj , (3. 14) 
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with real hg and Vq. This ansatz breaks the symmetry of spatial translations as well as the 
chiral rotations 



yq,n ^ y q^n 



3« en (?£ 



5,n; 



"'q,n "'q,in ^ "'q,n^ 



(3.15) 



and, as we discuss in Section III B[ implies a helical structure for the chiral condensate. 
A substitution of Eq. fl3.14p in Eq. (13. 6p gives 



hq = 2{d+l)aqVqJq{Q), 



where 



7,(Q) 



d + l 



1 + ^ cos (qQi) 



1=1 



and a subsequent substitution of that into Eq. (13.81) . with a use of 



dh 



4 X <^ 



nq 



2 
3 



3/1? „ + 2/i2,n q = 1 

2/ii,n g = 2 

q = 3 



gives the final form of the mean-field equations 



h2 
h3 



2 0373 

3a272 



2a272 

^373 



9ai7i / 2 , — z ,z ^ , , , 



hih 



3; 



0373 

16(ci + l)a373 
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X A{hg,Q;n). 



(3.16) 



(3.17) 



(3.18) 



(3.19) 
(3.20) 
(3.21) 



Also, ai,2,3 = i) ^) ""^j the coefficients in F3 (see Eq. (12. 5p ). and 71,2,3 are given in 



Eq. (13.171) (Here, for brevity, we write 7^ to denote 7g((5)). The function A{hq,Q;jj,) is 
defined as 

-3.Qne„ ^ C^tr log[S„l + D] 

and the derivative here is evaluated at 



A{hq,Q;n) = {hnhn) ■ e 



(3.22) 



= A[hl + 2hih2 + -h3 



(3.23) 
(3.24) 
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To determine which of the solutions of the mean-field equations is the true vacuum we 
calculate the mean-field free energy. Using Eq. fl3.16p . the energy per number of sites Ng is 
given by can be expressed in terms of the fields hq alone, and we find 

3/^2 

where the contribution of the b fields is 

^matter = "^tr log [S„l + D] . (3.26) 

We calculate A{hq,Q] jj) and £^matter in Appendix [XJ where we also show that the right 
hand side of Eq. fl3.22p is indeed independent of the lattice site index n. The result for 

^^matter IS ^ 

^matter = E / i^b - 3/i) 9 {E, - 3/i) , (3.27) 

where 9{x) is the step function. The momentum integrals are from — vr to vr, and E±, that 
replace £'(^)'(2) of Eq. fl2.3ip . are the band energies of the baryons. In this paper we focus 
on (i = 3 and Q = Qz, for which the dispersions are 

sinh^^i = sin2(p^/2) + sm\py/2) + (^■^J^Tsl± \c,\^^ , (3.28) 
s, = sin(p^/2)cos3g/2, (3.29) 
Cz = cos(p^/2)sin3g/2. (3.30) 

(For the general dispersions with Q = {Qx,Qy,Qz) see Appendix lAl) . S is related to the 
expectation values Vg = (ml) through Eqs fl3.6l) - fl3.8p and for Q = it is the spectral gap 
that defines the baryon mass ttib (see Eq. f l3.28p ) 

sinhm^ = E. (3.31) 

From here on we refer to energy bands of Eq. (13.280 with S > as massive bands, and to 
energy bands with S = as massless bands. 

In Fig. [2] we plot E± for p = (0, 0,^2). In the upper panel we choose = 3 (or S ~ 10), 
which is close to the value that solves the mean-field equations in c? = 3 and fi = (see next 
section). We present the cases of Q = 0, where E± are degenerate, and Q = n/S, where this 
degeneracy is removed. In the lower panel we show the case of = (or S = 0), where 
E± are again degenerate and also independent of Q. 
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FIG. 2: The energy bands E± of Eq. (j3.28p for p = (0,0, p^). Upper panel: Massive bands 
with rriB = 3 and Q = vr/S (highest and lowest curves), and Q = (middle curve). 
The bands E±{Q = 0) are degenerate and very flat, ranging between E = 3 and E = 
sinh^(3) + max ^X]f=i sm'^{pi/2)j ~ 3.0147. Lower panel: The massless band with S = 0. 



Finally, to make a connection with the discussion in Section IIIBI leading to Eq. (12.461) , 
we write Eq. (13.271) as 

i^matter = + ~ 3finB, 



where the anti-baryon and baryon energies, —\Sb\ and £b, are given by 



and the baryon number density is 

1 



(3.32) 



(3.33) 

(3.34) 



(3.35) 
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IV. THE MEAN-FIELD GROUND STATE 



In this section we present the solutions for the mean- field equations Eqs. fl3.19l) - fl3.21l) . 
and their corresponding energy densities. To understand the results we then study a formal 
limit of Eq. f l2.12p where one takes ^ 1. 



A. Physical case : N = 3 and d = 3 

We begin by describing the solution to the mean-field equation that leaves chiral symme- 
try intact in Section [I V A 11 and proceed to discuss the chiral broken phase in Section HV A2[ 
Finally, in Section IIV A 31 we discuss how the ground state evolves with increasing values of 



1. Solution I : Intact chiral symmetry 

The first solution to the mean-field equations is obtained by taking hi = h2 = = 0. 
This corresponds to a phase with intact chiral symmetry, Vi = V2 = V3 = 0, and zero baryon 
mass S = sinh(mB) = 0. In this case the energy bands are trivially independent of Q (since 

^mean-field _ degenerate 



E+ = sinh 



-1 



. ^sin^(p,,/2) = -Emassless- (4.1) 



The energy and baryon number densities of this solution are independent of Q as well, and 
are given by 

(dp \ 
j (-E'massless " 3/i) (-^massless " 3/i) , (4.2) 

Ub,! = J 6 (3/i - ^massless) • (4.3) 

In Fig. Owe plot Sj as a function of /i, where one can see that for /i = 0, we get Si ~ — 1. 
As we increase fi, the energy density increases, and saturates when /i = ^ 0.439, where 
it reaches zero. Looking at the form of -Emassiess from Fig. [21 we can understand why. The 



maximum energy to which -Emassiess of Eq. (14. ip can reach is sinh ^ YmaxX]?=i sin^(j9j/2) ~ 
1.317. This means that when fi > fig, where the baryon chemical potential 3fi is already 
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larger than 3 x 0.439 ^ 1.317, the band -Emassiess is aheady saturated. This saturation can 
be seen also in the plot of the baryon number density, ubj which reaches ubj = 1 at fig. 




FIG. 3: The energy density, £i, (solid line) and baryon number density, ubj, for the massless 
energy bands of solution type I, as a function of the chemical potential fx. 



2. Solution II : Spontaneously broken chiral symmetry 

In this case /ii,2,3 are nonzero, and we divide Eq. (13.191) with hi to obtain 

2 0373 1 20272, 

hi = hs. (4.4) 

9 ai7i % 0373 

Substituting this in Eq. fl3.20p and Eq. fl3.24p we get S as a function of h^, which we then 
substitute into Eq. (I3.2ip and obtain a single equation for h^, that we solve for all values of 
Q and /i. 

The results are presented in Fig. H] where we show niB = sinh^^(S) as a function of Q, 
for /U = 0, 0.5, 0.75, 0.9. In Fig. [3] we show the energy densities of these solutions. There are 
two branches in both figures that correspond to two solutions of type II that we find. From 
Figs. HHHl we see that at ;U = the minima of both branches (referred to as branches no. 1-2 
below) is at Q = and have 

mgVs ^ 1.0347, with Sff ~ -1.6826, (4.5) 
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FIG. 4: The values of ^msiQ', fJ-) that solve the mean-field equations. Note the two branches of 
the solution, and that one of them (branch no.l, ending at Q = 0, and ms/S = 1.0347) hardly 
changes with /i. We find that this branch has a lower energy and therefore corresponds to the 
ground state (see text). 



mgVs ^ 0.1309, with 



(2) 



-0.8981. 



(4.6) 



Consequently, the ground state of Solution II at yU = is given by the minimum of branch 
no. 1 (Eq. (S3D). 

To see whether the chiral wave instability takes place, we consider Fig. [5] and find that 
the minima of Sjj is at Q = for all values of /i. This is not in contradiction with the 
discussion of Section IIIBI because, there, we only showed that the baryon contribution to 
the free energy, £^b(S, Q), has a minimum at Q 7^ 0. What we see here is that the remaining 
contributions to S offset this minimum and lead to a ground state with Q = 0. 

It is important to note that even if the minimum of Sb was deeper, such that Su itself 
had a minimum at Q ^ 0, it is still unlikely that the chiral density wave would become the 
ground state. To understand why, recall that in Section [11 BI we showed that what allows the 
crystalline phase to compete with the S = phase of Solution I are the gapless excitations 
in the baryon spectrum. From Eq. (13.281) we see that the momentum of these excitations is 
given by 



P = (0,0,pj, 
19 
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FIG. 5: The energy density Eji as a function of Q and ;U for solution II. The lower branch (branch 
no.l) for each fi is the one that ends at Q = and m^/S = 1.0347 (see Fig. H]). 



which has a solution only if 



sin2p^/2 = sin2(ArQ/2) - 



E < 1. 



(4. 



(4.9) 



Unfortunately, with the baryon mass of our ground state (14.51) we have E = sinhm^ ~ 11. 
This prohibits gapless excitations, and makes the appearance of the chiral waves improba- 
ble. ^ 



3. Evolution of mean-field ground state as a function of 

Here we discuss the way the ground state evolves with increasing chemical potential /i. 
Given the conclusions of the previous sections, we restrict to states with Q = Q only. 

At yU = 0, solution of type I has £i —1, while the solution of type II that has the 
lowest energy is given in Eq. (I4.5p . Comparing Sj and the energy in Eq. (14. 5p we see that 
the latter is the ground state of the system, and so, the baryon mass is m^/S ~ 1.0347. The 

The true relevant value of S that should tested against Eq. (|4.9p is S(Q;/i) = sinhmsiQ', n), but from 
Fig. H we see that S((5; fj,) > E(7r; 0) ~ 4, for which Eq. (|49|l is stih not satisfied. 
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condensates Vq (Eq. (12.181) ) that we find for this ground state are 

Vi ~ 0.67685, 



) 



XXY 
3 J 



Vo 



0.31947, 
0.07872. 



(4.10) 
(4.11) 
(4.12) 



Comparing these results to the hterature, we see that our mass is 0.22% — 0.54% lower than 



other tree level results (see, for example, the summary in 
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23 
21 



), and that our condensates 



221 . |23| (discrepancies are on 



are close to the analytical and numerical results of 
the level of 2% -5%). 

When n increases, the energy densities of both type of solutions grow. As explained 
in Section IIV A 11 the energy of solution I stops increasing at where it reaches Sj = 0. 
Nothing special happens to solution II at this point, and its energy is still negative and 
continues to grow monotonically. In Fig. [H] we replot Sj{fi) of Fig. [S] together with the 
ground state value of Sjj{fi) (i.e. for given /x, we extract the minima of the curves in Fig. E]). 




FIG. 6: The energy densities £i and En as a function of /x. 



From Fig. [6] we see that a transition occurs at fi = fit where 



Siiif^t) = 0. 



(4.13) 
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SubstitutingEq. and Eq. (IHTTUD in Eq. flTO|) gives (for Q = 0) 

^{d+l)j2a,V; = j - 3/i,) ^ - . (4.14) 

From Fig. Hlwe see that, for (5 = 0, the mass does not change as a function of /i up to 
yU ~ m^/S, at which point solution II disappears. This reflects the fact that the condensates 
Vq are independent of /i as well, and allows us to use the /i = values of and Vq in 
Eq. fl4.14p . This is a consistent procedure provided that the resulting ^it obeys yUt < rriBl'i- 
Indeed, solving Eq. (14.141) for /ij gives 

/it = {^^\ (i)'^(P)) - ^[^iVl ^ o.,Vl ^ a,vi) 

~ 1.0367- (0.45813 + 0.01914- 0.001452) = 0.5609, (4.15) 

which is lower than itlb/^ — 1.0347. Consequently this means that the chiral symmetry 
of Eq. (II. 3p is restored in a first order transition at — 0.5609 that separates a low 
density phase characterised by the baryon mass of Eq. (14. 5 p and the chiral condensates in 
Eqs. (I4.10p - (l4.12p . and a high density phase where this symmetry is intact and the baryons 
are massless. 

Let us emphasise that the way chiral symmetry is restored here is largely influenced by 
lattice artifacts of the infinite coupling limit. In particular, we cannot observe the crystalline 
phase partly because the mass of the baryons in lattice units, m^, is around ~ 3. This 
prevents the baryons spectrum from having gapless excitations, and makes their contribution 
to the free energy SBifriB, Q) higher compared to what it would be with = Q = 0. This 
is in contrast to the 'continuum- like' scenario that we discussed in Section [11 B[ where 
was small, and EsijnBjQ) < £^^(0,0). In the next section we study such a continuum- like 
case by taking a formal limit of ^Shks; iii which the mass of the h fields is significantly smaller 
than discussed above. 

B. The formal limit of iV ^ 1 

In this section we study a formal limit of Eq. (I2.12p . where one takes = 1 and Fn{u) = 
u/4. In this case Eq. (I2.12p cannot be obtained from an underlying lattice gauge theory, 
but nevertheless seems to be useful to understand, since it shows chiral restoration with less 
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lattice artifacts. In their original preprint, the authors of [15| have considered this limit 
as an example for a case where the baryonic terms in Eq. (12.121] are important at /i = 0. 
Here we approach this limit at nonzero fi and seek for signs of a crystalline solution for the 
mean-field equations. We first do so for c? = 1, where it is known that crystalline instabilities 
are robust [2, l3[. (The physical significance of the result in this case is, however, unclear 
since the continuous U{1)^ symmetry cannot break in 1 + 1 dimensions.) 

We begin the analysis by noting that for = 1, Eq. (12.121) depends only linearly on 
m„m„+p. This allows us to take h2 = = V2 = = 0, and to call hi = h,Vi = V. The 
mean-field equations become (we keep d general at this stage) 

h = 2{d+l)aiV-fi{Q), (4.16) 
V = A{h,Q;fx), (4.17) 

and the free energy is now given by 

4(d+ l)ai7i((5) 

Both A{h,Q;fi) and £^matter are still defined in Eqs. (I3.22p . and (I3.26p . with the difference 
that the dependence of S„ on /i„ is now 

= hn = he-''"^''K (4.19) 

Again, it is easy to see that there are two type of solutions which are the analogues of 
the type I and II solutions in the SU{3) case. In Fig. [7] we present the energy bands for 
niB = sinh~^(S) = 0.5, which is close to the value that solves the mean-field equations for 

= 1 and d = 1. With this value of the effect of the lattice coarseness is smaller 
compared to the A^ = 3 case : the massive and massless bands of solutions I and II overlap, 
and there exists momenta where the energy of solution II is zero. 

As a result, a crystalline ground state is stable. In Fig. [H]we present the dependence of 
V and Q on /i. We also plot the following analytic expectation for Q{^) that arises from 
our discussion in Section Hi Bl (adjusted to A^ = 1) 

g(^) = pp{^) = 2 sin"^ (sinh /i) . (4.20) 

It is clear that Eq. (14.201) works well, especially for S ^ 1, which means that the mechanism 
described in Section III Bl is indeed the one generating the chiral density wave instability 
here. 
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FIG. 7: The energy bands E± of the baryons for p = (0, 0,^2)- The dashed hne is the gapless band 
with = 0, and the sohd hnes are have ms = 0.5, and Q = 0, vr. 

When we move from d = 1 to d = 2 the crystalline phase disappears. We find that the 
reason is that the meson self energy described by the first term in Eq. (14.181) prefers Q = 0, 
and 'wins' the instability generated in £^matter- This effect comes from the function 71 (Q) in 
Eq. (13.171) . and can be considered as a lattice artifact as well, since in the continuum limit the 
lattice quantities and Q vanish, and so 71 (Q) 1. Indeed, in continuum treatments 



of the Nambu-Jona-Lasinio model such as 



2J], the contribution to the free energy that is 



the analog of the first term in Eq. (I4.18p . does not depend on Q. 

When we move to = 3 it appears that even without the Q dependence in 71 (Q), 
the crystalline phase is unstable. This changes if one makes the meson self-interaction 
Fi{u) = \u stronger by changing it to Fi{u) = aiu with ai > 1/4. 



24 



o Z 




0.2 0.4 0.6 0.8 



FIG. 8: The values of V (circles), and Q/ir (dots) that solve the mean-field equations (|4.16|) - (j4.17|) 
in 1 + 1 dimensions. The red solid line denotes the analytic expression of Eq. ()4.20p . Note that we 
present Q/tt (black square symbols) only when S > 0, since otherwise gj^'^^^-^'^^'^ = g, and Q is 
ill-defined. 

V. CONCLUSIONS 

Our aim in this paper was to see whether the phase diagram of strongly-coupled lattice 
QCD at zero temperature and nonzero density includes an incommensurate crystalline phase 
with chiral density waves. To do so we followed Hoek, Kawamoto, and Smit (HKS), and used 
the effective hadronic action, ^rks, that they derived in [15] for one-component staggered 
fermions. We formulated a mean-field theory for Shks which is novel in the following ways. 
Firstly, we do not neglect the terms in ^kks that have four and six powers of the meson 
fields, and treat the baryon contribution in full. Secondly, our mean-field analysis is free of 
any undetermined parameters and we fix all vacuum expectation values by minimising the 
free-energy. Lastly, we do not assume a homogeneous ansatz, which is imperative in order 
to look for the crystalline phase. This leads us to introduce auxiliary fields without the 
usual Hubbard-Stratonovich transformation, that can become ill-defined for inhomogeneous 
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vacua. 

Despite the fact that the generic structure of ^kks can give rise to a crystal phase (see 
discussion in Section [II Bl) . we find that this does not happen with our mean-field ansatz. 
This is partly due to lattice artifacts that come in the form of a large lattice baryon mass 
~ 3. The latter prohibits any possibility to have a gapless baryonic spectrum in the 
presence of the chiral density waves, which would allow the crystalline vacuum to compete 
with the massless homogeneous vacuum. To check this explanation, we take a formal limit 
of S'hkS) where itlb — 0.5, and indeed find that this lower mass can stabilise the crystal. 

More precisely, when we study this formal limit in the 1 + 1 system, we see a crystal 
structure at intermediate densities, with a wave vector given by the Fermi momentum of 
the system. However, this result is not robust; when we move to 2 + 1 and 3 + 1 dimensions 
we see that the crystalline instability is too weak to survive. The reasons for that include 
an additional lattice artifact, that does not go away when decreases, and a too weak 
interaction between the fermions. 

Our results stress that to get continuum physics from Monte-Carlo simulations at nonzero 

and low temperature, one will have to simulate on relatively fine lattices, with weak 
couplings. In particular, the baryon mass in lattice units, m^, should be much smaller 
than its strong-coupling limit of ~ 3. Only when this happens will the structure of 
the baryon energy bands be 'continuum-like'. For stronger couplings, the energy bands are 
qualitatively different than the continuum ones, and physical phenomena that are sensitive to 
their structure will also be qualitatively different than it is in the continuum. Chiral crystals 
is an good example for such a phenomenon, that simply goes away at strong couplings. 

There are several ways in which one could extend this study. First, recall that we have 
focused on an ansatz that breaks translation invariance only in one direction. This direction 
was chosen to be along one of the lattice axes, but a further study can generalise this easily 
(we have already derived the free energy for a general Q = (Qi, Q2, Qs)) and check whether 
this makes the crystalline phase more robust. Second, as the authors in [4i] point out, 
one can lower the energy of the crystalline phase by considering an ansatz which is a linear 
combination of waves, that will break translation invariance in all directions. Third, it will be 
useful (but also hard) to see how do the sub-leading terms in the strong-coupling expansion 



influence the results. These were obtained in 



25l | for the fermions we consider here, and give 



rise to more complicated interactions between the hadrons. Finally, it can be interesting to 
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numerically study the possibility of cliiral density waves in the Nambu-Jona-Lasinio model, 
where Monte-Carlo simulations are free from the sign-problem. 
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APPENDIX A: THE BARYONIC DETERMINANT 

In this section we calculate the baryonic determinant Eq. (13.261) and its variation 
A{hq,Q]fi) of Eq. fl3.22p . We begin by showing that A, which, for the general SU{N) 
case, is defined by 

A{h„Q;f^) ^ (6„60e-^*«"^", (Al) 
- / DbPb Kbn 

^^"^"^ = JDlDb eso 

So = J2 ^" + ^rim] bm, (A3) 

n,m 

is independent of the lattice site index n. To see this we use the symmetries of the action 

So = bKb = Y.K Se-^^Q"^" bn + W'^n {Vnubn+u " VnX^u) = Sj: + ^kinetic, (A4) 
71 nu 

where = ± is the parity of the site. In particular, S'kinetic is invariant under spatial 
translations 

: n^n + R (A5) 

and chiral symmetry 

f/(l), : 6„^e^^^^"6„, 6„->6„e*^^^". (A6) 

The term 5*2 is invariant under these symmetries only for the correlated transformation with 
26 = Q ■ R. (More precisely, if R takes one from an even site to an odd site (e„ = —en+R) then 
an extra transformation of reflection around the origin is needed to keep Sj: unchanged.) 
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Although 5*0 is invariant under these correlated transformations, {bnbn) is not, and trans- 
forms as follows 



For er. 



For 



T- 



reflection iv 



{bnbn) K+Rbn+Rl K+Rbn+li) 



-NiQRen 



(AT) 



rbn) ^ {K+nbn+R) (K+R bn+ll) ^ 

-NiQRe„ _ 



NiQR€„ 
+NiQRe_ 



Putting n = gives = (^_r^k) e'^^'^^^R for = +1 from Eq. (]A7[) . and for = — 1 

from Eq. ( lASp . This proves that 

(bnbn) e-^^'^n.n _ independent of n, (A9) 

and therefore that^ 



A{hg,Q;fi) 



1 d log det [S e"^*^"^" 1 + D 



We now turn to calculate the determinant itself. It will be convenient to note that 



(AlO) 



log det K = log det 



and to calculate det K = det 



log det [1 + S-^e^^^"^"D]+^logS-Arz^Qf?e„, 

n n 

(All) 

1 + S-ie^*'5"''"D instead. (Note that the last term in 
Eq. (lAlip drops out for a lattice with an even number of sites in each direction). To proceed, 
we use coordinates defined on a new lattice with spacing a = 2, and write n = 2X + p, with 
Pi,...,d = or 1 and X taking values in the new lattice [l9|]. In these coordinates we use the 
following definitions 



bn = bp{X), 

bn+v = J2(^p',p+^bp'{X) + 5p'^p^0bp>{X + u)) , 
p' 

bn~o = ^{5p',p^0bp>{X) + 5p'^p+obp>{X - u)) . 



Moving to momentum space with 



bp{X) 



2^ 

—y 

Yd 



MpX+epQi2X+p)/2 



bp{p), 



bp{X) = J^ye-<^''^^^'^^''^^'^/%{p), 

V -'*s p 



(A12) 
(A13) 

(A14) 

(A15) 
(A16) 



^ This is different from Eq. p.22p since here we evaluate the determinant for the hehcal ansatz, and then 
take the derivative with respect to S. 
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we get 



^0 = bKb = j2b,ip)k^,ip%{p), 



(A17) 



p' = p- NQe 



pi 



(A19) 



where we used ep = — Cp/. Eq. ( lAlSI) can be brought to the form 



K{p) = 1 + zS-i E sin (p,/2 - iNij6,^o + NQ,e/2) (A20) 

= 1 + zS"^ [ro(cj) sin {uj/2 - zA^/i) 

+ E r.(p) sin(p,/2) cos(ArQ,/2) + T,{p) e cos(p,/2) sin(ArQ,/2) 



Here we defined 



and also 



where fji, = rjy for z/ G [1, (i] and ?7o = 1- The matrices Vy obey 



with 



that also obey 



{r,, r J = 25,A, {r., e} = {r,, Ts.} = o, 



(r5.)pp, = ^UP) iSp,p'+, - 5p,p'-.) e^^(''-''')/^ 



(A21) 



(A22) 



(A23) 



(A24) 



(A25) 



{Is., Ts.} = 25,^1, {Is., e} = 0. (A26) 
Next we use the fact that det e K{p) e K{p) = det K'^ip) and that {Vy{j)), e} = to get 

(det k{p)f = det ek{p) e x det k{p) = det (p), (A27) 



with 



K2{p) = {ekip)e)k{p) 



1 1 + s-M E - E + E SuC.T^r, e 



\u=0 u=l 



(A28) 
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Here we defined 



So = sin (po/2 - iN^x) , 



(A29) 



and for 1/ > 



= sin(p^/2)cos(A^(5^/2), 
= cos(p^/2)sin(Arg^/2). 



(A30) 
(A31) 



To proceed we need to find a matrix that anti-commutes with Fj^F^e for all i/ ^ /i. We can 
choose any of the matrices M = T^^ (which also obey = 1) and we have 

,4 



(det K{p)) = (det K2{p)f = det MK2{p)M x det K2{p) = det K^{p), (A32) 



with 



K^{p) = {MK2M)K = 



Using the anti-commutation relations of Fj, we get = —B ■ 1 with 



(A33) 



(A34) 



and defining D{p) = detir4(p) we obtain the following form 

D{p) ■ = (sin2(po/2 - iNij) + (p)) x (sin2(po/2 - iNij) + ei(p)) . (A35) 

Here 

el = |sl2sin2 0+ (^^E2 + |s12cos2 0± Icl^ , (A36) 

Si = sin(pi/2) cos(Arg./2), (A37) 

a = cos(p,/2)sin(Arg./2), (A38) 

COS0 = c • s. (A39) 

Since det K{—p) — det then det K = Y\p det ^(p) is reaP and we can write det K — 



^ This can be proved as follows. Wc have K{p,Q) = K*{—p,—Q), but also det K{p,Q) 
detT5i,eK{p,Q)eT^f, = det K{p,-Q). This means that det K{p,Q) det K{-p,Q) 

det K{p,Q) det K*{p,-Q) = det K{p,Q) det K*{p,Q) = \ det K{p,Q)\^ is real. 
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(det K'^y^^ = (UpdetK^ip))'^'' or 



1/4 



d+1 



logdeti^: = \T.^og{D{p)f = ^ £^ (^gj logD(p). 



(A40) 



Dropping irrelevant constants, we get 



-^logdeti^ = 7 y 



27r 



log sin2(po/2 - iNy) + 



1 91ogdeti\: E 



E 



f dp\ dpo 



2tx 



1 + 



k1 



2 _|_ 5^ ^-->c2 - 



^ 27r sin2(po/2 - iN ii) + e^(p) ' 



(A41) 
(A42) 



Performing the change of variables z = e*^" the integral in Eq. flA42p becomes a 
simple contour integral and using complex analysis one can show that 

r 1 ^ 0{E,{p) - AT/i) 

where sinh^ -^^(p, s) = e^(p). The integration in Eq. (lA4ip is carried out as follows (dropping 
an irrelevant constant of integration) 

1 



— log sin^(po/2 - iNfx) + el{p) 
-IT zvr L 



Xo J-7T 2n sm'^{po/2 - iNfi) + X 
2 [Ehip) - Nfi] 9{Eh{p) - N^i). (A44) 



The final result is then 

1 (91ogdetA: _ S 



E 



, , .d (i + b , 

2^/ ^ 



9iE,ip)-Nfx) 



-^logdetir = - y 

2£i 



dp' 
2^ 



{E,{p)-Nfi)9{E,{p)-Nfi). 



(A45) 



(A46) 



which we evaluate numerically. In particular, in the case of d = 3, we present results for 
Q\\z, in which case 

|2 



ei = |s^r+u/S2 + sl±|c| 



s, = sin(p,/2)cos(A^g/2), 
c = cos(p2/2)sin(Arg/2), 
= sin^pa,/2 + sin^Py/2. 



(A47) 
(A48) 
(A49) 
(A50) 
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Since the transverse directions appear very simply here, we define the density of states 

D(s)ds = — ^ sin^(pj,/2) — s for the transverse momenta. One can show that 

D{s) = ■^K{s{2 — s)) where K{x) is the complete elliptic integral of the first kind (for 
example see |26|). 
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